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Abstract 

Stationary states of non-interacting electrons on a Koch fractal are investigated within a tight 
binding approach. It is observed that if a hierarchically long range hopping is allowed, a suitable 
correlation between the parameters defining the Hamiltonian leads to spectacular changes in the 
transport properties of finite, but arbitrarily large fractals. Topologically identical structures, that 
are found to support the same distribution of the amplitudes of eigenstates, are conducting in some 
cases and insulating in the others, depending on the choice of the hierarchy parameter. The values 
of the hierarchical parameter themselves display a self-similar, fractal character. 

PACS numbers: 73.21.-b,73.22.Dj,73.23.Ad. 
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I. INTRODUCTION 

The localization of single particle states in a disordered 
system, viz. the Anderson localization has been ex- 
plored extensively since its inception way back in 1958. 
The investigation didn't remain confined to the cases of 
randomly disordered systems in one, two or three dimen- 
sions alone, but also encompassed quasi one dimensional 
systems with topological disorder t^J-iSL quasi-periodic 
chains and deterministic fractals |9|- jl8|. The ba- 

sic result is that, all the single particle eigenstates of one 
and two dimensional randomly disordered lattices should 
be exponentially localized irrespective of the strength 
of disorder, while in three dimensions there is a metal- 
insulator transition. 

Exceptional results however, started surfacing in the 
last decade of the twentieth century, when it was shown 
by Dunlap et al that a certain kind of correlated dis- 
order could lead to a finite number of resonant eigenstates 
in an otherwise disordered one dimensional (1-d) chain. 
Such states are extended but, of non-Bloch character, 
and lead to a ballistic transport of electrons at special 
values of the energy. The so called random dimer model 
(RDM) proposed by Dunlap et al [l^ led to a thorough 
investigation of the effect of correlated disorder on the 
electronic and other excitonic states in one dimensional 
systems [20| - 27 1, where, ordinarily an un-correlated dis- 
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order renders all the states exponentially localized. 

The effect of such short range correlations among the 
constituents in a one dimensional system has also led to 
exotic distribution of extended eigenstates in one dimen- 
sional quasi-periodic lattices [28j-|3T| where, owing to the 
self similarity of such systems, one encounters an infinity 
of such extended states in an infinite 1-d system. 

Deterministic fractals have also played their part. A 
subset of these structures has been shown to possess ex- 
tended non-Bloch single particle states with a variety of 
transport properties |32j-(37|. In contrast to the RDM 
or the quasi-periodic chains in 1-d, the fractal lattices do 
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FIG. 1: (a) A 2nd generation Koch fractal with long range 
hopping. The sites B{j) with j = 1 and 2 denote the 'bulk 
sites' where j denotes the stage of hierarchy. Accordingly, 
the longer ranged hopping integrals are assigned values T{j), 
j = 1,2. The boundary sites are marked as A and the re- 
maining sites are C. (b) A unit cell (inside the dashed box) 
of a periodic approximant built by repeating a first generation 
Koch curve is shown, is enclosed in the dashed box, and, (c) 
The final trimer that is used for the trace analysis. 



not exhibit any short range clustering of atomic sites (in 
the sense of the RDM) that could possibly lead to lo- 
cal resonances resulting in a de-localization of the wave 
function. Rather, the topology of the entire lattice plays 
the key role in giving rise to such extended eigenstates 
at special values of energy of an electron hopping around 
in these lattices. 

The extended states in the cases referred to, appear, in 
general, in isolation in the energy spectrum where they 
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coexist with an otherwise singular continuous part. How- 
ever, in recent times extensive numerical studies have 
pointed out towards a possibility of having a continuum 
of such extended states in a class of deterministic frac- 
tals 13- [13. 

A universal feature of such extended states, whether 
they arise in RDM models or in the fractals reported so 
far is that, their existence is not sensitive to the numer- 
ical values of the Hamiltonian parameters. The short 
range clustering of the atomic sites or the overall lat- 
tice topology determined the criteria of the existence of 
unscattered eigenstates. In the present communication 
we report an exceptional case where, well defined corre- 
lations between the values of a subset of the Hamilto- 
nian parameters is essential for the construction of an 
extended eigenstate which has equal amplitude at every 
lattice point. We explain the phenomenon by solving 
the Schrodinger equation on a Koch curve with hierar- 
chical interaction [38]. A Koch curve with a long range 
hierarchical distribution of electron-hopping models, in 
an approximate manner, a polymeric system [ssj . We 
show that this model can lead to unusual spectral and 
electronic transport properties when the hierarchy pa- 
rameter is allowed to vary, mimicking (in a very crude 
sense though) local distortions in the system. The most 
astonishing result is that the same distribution of eigen- 
states, which makes the lattice conducting for a given set 
of the Hamiltonian parameters, represents a completely 
localized state with zero transmission across the Koch 
fractal. In addition to this, we find that such contrasting 
results are obtained for an infinite variety of values of the 
hierarchy parameter. These values are distributed in a 
self-similar manner, as is apparent from a study of the 
two-terminal transport across a Koch fractal. We work 
within a tight binding formalism and employ a real space 
renormalization group (RSRG) decimation scheme to un- 
ravel these strange states at all scales of length, and to 
obtain the end-to-end transmission properties of a Koch 
fractal of any size. 

In what follows, we explain our work. In section II the 
model and the method are described. Section HI contains 
the results and the discussion, and in section IV we draw 
our conclusion. 



FIG. 2: Eigenvalue spectrum of a Koch fractal as a function 
of the hierarchical parameter A obtained from the trace of the 
transfer matrix for a 7th generation fractal, taken as the 'unit 
cell'. 



illustrate our results for the specific case when the longer 
ranged hopping integrals assume values T{n) = R'^t, R 
being the hierarchy parameter and n represents the level 
of hierarchy. For a Koch fractal of finite but arbitrarily 
large size, the on-site potential will be assigned val- 
ues sa (the extreme sites) and 63(72), £c(n) etc for the 
bulk sites depending on their positions in the lattice as 
explained in Fig. [T] 

Owing to the self similarity of the lattice an applica- 
tion of the RSRG decimation method is used to scale 
the original system. An appropriate subset of sites are 
decimated leading to the recursion relations of the Hamil- 
tonian parameters, viz. 
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II. THE MODEL AND THE METHOD 



where, 



We start by referring to Fig. [TJ 

The tight binding Hamiltonian in the Wannier basis 
reads. 
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tiic]c^ + h.c. 
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where, Ci is the on-site energy of an electron at the site 
i. tij is the hopping integral between neighboring sites. 
For nearest neighbor sites along the ridge of the Koch 
curve we assume tij = t, while the longer range hopping 
integrals can, in principle, take on any values. We shall 
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The above set of recursion relations will be used to ex- 
tract information about the eigenvalue spectrum and to 
discern the nature of the eigenstates in the light of the 
motivation presented before. 
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FIG. 3: (a) Renormalized hopping between the extremities of 
a 9th generation Koch fractal against the hierarchy parameter 
R at an energy E = eA + t- (b) A fine scan of the values of R 
reveals the fractal nature of the distribution of the values of 
R. 



by sequentially applying the recursion relations Eq. [2] n- 
times, bring it down to an effective first generation block 
and subsequently to an effective triatomic molecule as 
discussed above to get the trace of the transfer matrix. 
The spectrum of an infinite lattice is obtained when the 
unit cell itself becomes infinitely large. 

We present the result in Fig. [2] where the unit cell 
is a seventh generation Koch fractal with hierarchically 
distributed values of the long range hopping integrals 
T(n) = R^t. The 'band structure' is displayed as the 
hierarchy parameter R is varied. The behavior seems to 
saturate, at least becomes indistinguishable, for unit cells 
of larger size. So, the present figure can be taken to rep- 
resent a very good approximation to the true eigenvalue 
spectrum in the infinite lattice limit. 

With low absolute values of the hierarchy parameter 
i?, the spectrum exhibits clustering of eigenvalues, giv- 
ing rise to the formation of sub-bands. The sub-bands are 
separated by gaps which widen as R increases. However, 
for R quite large compared to the nearest neighbor hop- 
ping t, the sub-bands shrink in their widths and global 
gaps open up. An indefinite increase in the value of R 
finds the spectrum consisting of a sparse distribution of 
points. 



III. RESULTS AND DISCUSSION 

A. The eigenvalue spectrum 

We try to gain some insight into the spectrum of an 
infinite Koch fractal in the presence of hierarchically dis- 
tributed long range interactions. To achieve this, we pro- 
ceed in the spirit of the one dimensional quasi-periodic 
lattices The basic idea is to repeat a finite gener- 

ation Koch curve periodically, pick up the unit cell and 
extract the energy values for which the trace of the trans- 
fer matrix for the unit cell remains bounded by two Q . A 
prototype unit cell for a periodic approximant generated 
by repeating the first generation Koch curve is shown in 
Fig. [IJb) . To construct the transfer matrix and work out 
its trace in this case, the top vertex C is decimated to 
reduce the system into a 'triatomic molecule' A1-A2-A3 
(Fig. (He)) where, 

(■Ai — 



and, eA-i — £^2- The 'inter-atomic' hopping integrals are 
tA,-A2 = t, and tA,-A, = T(l) + ty[E - ec{l)]. The 
transfer matrix across this triatomic molecule is then eas- 
ily calculated. One can, in principle, begin with a fractal 
at any arbitrary generation, renormalize the parameters 




FIG. 4: (a) The transmission coefficient of a 7th generation 
Koch fractal with hierarchical interactions — 1 < -R < —0.6, 
and (b) a finer scan for —0.83 < R < —0.71. The fractal 
distribution of T against R is to be noted. 
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FIG. 5: Transmission coefficient of a 6th generation Koch 
fractal against the hierarchy parameter R. Here, the energy 
has been obtained from the equation E = e'^ + t' , that is, 
from a one step renormalized lattice. The variation with R 
and its fine scan revealing again a self similar character, are 
shown in figures (a) and (b) respectively. 



B. The atypically extended eigenstates 

In this discussion we highlight certain aspects of a sub- 
set of eigenstates on a Koch fractal of finite but arbitrar- 
ily large generation. Let us choose the energy of the elec- 
tron (the Fermi energy) as E = eA+t. With this choice, 
it is easy to verify that one can construct an eigenstate 
whose amplitude at each site of the Koch fractal will be 
identical (normalized to unity, say) if we select 



ec{n) = €A-t 
esin) + T{n) ^ SA-t 



(5) 



for all n > 1. The above choices of the on-site poten- 
tials and the hierarchical hopping integrals T{n) define a 
correlation between the values of the parameters of the 
Hamiltonian. Interestingly, we can see that once the on- 
site potential ca of the extreme sites and the nearest 
neighbor hopping matrix element t are fixed, €3(71) and 
T{n) can, in principle, be chosen from a set of random 
values, but in such a way, so as to satisfy the desired 
correlation, viz, e_B(n) + T{n) — sa — t. Since, with this 
choice, the probability of finding the electron is same 
on every lattice point, it is not unusual to assume these 
states to be extended in the conventional language. How- 
ever, a study of the flow of the nearest neighbor hopping 
integral t under successive RSRG steps reveals a deeper 
story. It is known that, if, corresponding to certain value 



of the energy E, the hopping integral remains non-zero 
under RSRG iteration, it implies that the nearest neigh- 
bors are 'connected' at all scales of length. Hence, the 
corresponding eigenstate can legitimately be designated 
as an extended state. In the present case, the result of 
this investigation is found to be strongly sensitive to the 
initial choice of the parameters e_B(ri) and T{n). 

To clarify the above idea, we consider the particular 
situation where the longer range hopping integrals T{n) 
are distributed in a hierarchical fashion, viz, T{n) = R^t. 
We set E = ca + t, and take €c{n) — €a — t for all 
n. The choice of £3(71) is then restricted by the desired 
correlation as it appears in the last of the equations in 
Eq. [5j A Koch fractal, on successive renormalization by 
an appropriate number of steps may be reduced to an 
effective diatomic molecule comprising of the end sites 
only. The effective hopping between the surviving end 
sites will decide on the extendedness of the wave func- 
tion at that particular energy. Most interestingly we flnd 
that, the flow of the hopping integral t under successive 
iterations is strongly sensitive to the choice of the hierar- 
chy parameter R (which automatically fixes £5(71.)). We 
present one such result in Fig. [31 where the end-to-end 
renormalized hopping in a 9th generation Koch fractal 
is plotted as a function of the hierarchy parameter R. 
We have taken E = I and T{n) = R"-t with = 0, 
t = 1, ecin) = —1 (for all n), and eB{n) = eA — t — T{n). 
It must be appreciated that for every value of i?, with 
the above conditions being satisfied, the wave function 
displays an identical pattern of amplitudes on all sites, 
yet, we find that the end-to-end hopping becomes zero or 
remains finite depending on the values of R (and hence 
the values of tB(n)). In addition to that, a finer scan of 
any of the intervals of R reveals that the values of the 
hierarchy parameter for which the end-to-end hopping 
across a Koch fractal of any generation remains non-zero 
(a true extended state), are distributed in a self-similar, 
fractal pattern. Every value of i?, defines a new point in 
the parameter space [e^, ^sin), ec{n), T{n), t]. therefore, 
for E = tA -'r t, we come across a fractal distribution of 
points in the parameter space. The features described 
above are reflected in the end-to-end transmission across 
a flnite Koch curve with hierarchical interactions. The 
results are presented below. 



C. The transmission coefficient 

We have examined the two-terminal transport of a 
Koch fractal of an arbitrary generation by placing the 
lattice between two semi-infinite perfectly periodic leads 
characterized by a constant on-site potential eg and near- 
est neighbor hopping integral t^. A fc-th generation frac- 
tal is reduced, by applying the decimation recursion re- 
lations Eq. [5] k times, to an effective diatomic molecule 
comprising of the extreme sites A — A with renormalized 

(k) 

on-site potential e\ and renormalized 'inter atomic' hop- 
ping integral t'^^^ . The transmission coefficient r is then 



5 



obtained from the well known formula [42 



4 sin qa 

(6) 



[Mi2 - M21 + (Mil - M22) cosqaY + [Afu + M22Y sm^ 



Here, 
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-M12 
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The lattice spacing a in the leads is taken to be unity, and 
q = cos~^[(£' — eo)/(2to)]- Throughout the calculations 
we take eo and to = 1- 

Fig. U (a) displays the transmission across a 7th gen- 
eration fractal ioT E = eA + t for a certain choice of the 
range of the hierarchy parameter R. For every value of R 
in this range, we have taken £3(71) = eA — t — T{n) which 
makes the amplitudes of the wave function same at every 
site (normalized to one). The spectrum consists of gaps 
and resonant transmissions indicating that, for certain 
values of R the 7th generation Koch network may become 



completely transparent to an incoming electron with the 
above energy, while that same electron finds the network 
totally opaque for other values of R. A fine scan of a 
sub-interval of R reveals that the values of R presenting 
the above feature are indeed distributed in a self-similar, 
fractal manner. This is illustrated in Fig. SJ^b) and is in 
agreement with the initial observations made based on 
the recursive behavior of the nearest neighbor hopping. 

One can construct eigenstates for which the amplitudes 
will be same (unity) at all the sites of a renormalized 
Koch fractal as well. The procedure outlined above for 
the lattice at the bare length scale is to be implemented 
on a renormalized version of the system. For example, 
we can set iJ = -I- 1' on a one step renormalized lattice. 
The eigenvalues of the desired states are then extracted 
by solving the polynomial equation 



S-F(e^,i,eB(l),ec(l),T(l)) -0 



(8) 



where. 



F = 



eaE^ + E [e 



2e 



eA (£b(1) + ecjl) + T(l))] + sa [£B(l)ec(l) 

2 - [eBil) + ec(l) + T(l)] E + £5(1)60(1) - 2t^ + ec(l)T(l) 



E^ - [ 



6c(l)r(l)]-t^6c(l) 



(9) 



The potentials and the hopping integrals appearing as 
the arguments of F can be chosen independently. It is 
to be appreciated that the roots of the equations are free 
from esin) and ec{n) and T{n) for n > 2. We shall 
now demand 6(7(1)' = — t' which leads to an equation 



J 



defining ec(2) in the original lattice, viz 
6c(2) = [6:4-<'-A]^ 



where. 



(10) 



[E-eciW 



E + Til)-eB{l) £;2-£;[ec(l)+T(l)] + 6c(l)[6B(l) + r(l)]-2t2 



(11) 



The energy E is chosen from the solutions of Eq. [51 This 
settles the initial value of 6c (2) we should begin with. 
To mimic the earlier situation in the original scheme (the 
bare length scale) we now set ec{n) = 6^(2), as obtained 



above, for all n > 2 in the original Koch geometry.. 

In addition to this, we should set 63(71)' +T(n)' = 6^ — 
t' for n > 1, which leads to the selection of the un-primed 
values of 6b(?t-) and T(n) such that, eBi'ri) + T{n) = (6^ — 



6 



t' — A)E=c'^+t' = ec(2), for n > 2 on the original (bare 
length scale) lattice. The entire scheme now works just as 
it worked in the original scale of length. The process can, 
in principle, go on and every time we extract the energy 
eigenvalue by solving the polynomial equation E — e^'' — 
fik) _ g fj-oni the k—th stage of renormalization, we need 
to pre-define a bigger parameter space to begin with. The 
energy of course remains independent of eB{n), ec{n) or 
T{n) for n > k. The scheme works, the correlations 
become much more complicated though. 

IV. CONCLUSION 

In conclusion, we have unraveled a set of atypical ex- 
tended states in Koch fractals possessing hierarchically 
long range interactions. The states have identical am- 
plitudes at all atomic sites in the lattice for mutually 
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